Enhancement of spin stiffness with dilution in a ferromagnetic Kondo lattice model 
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Carrier-induced ferromagnetism is investigated in a diluted ferromagnetic Kondo lattice model 
for several ordered impurity arrangements on square and cubic lattices, allowing for quantitative 
comparison with different theoretical pictures. The spin stiffness is found to be optimized with 
respect to hole doping concentration, exchange interaction strength, as well as spin dilution due 
to a competition involving magnitude of carrier spin polarization and its oscillation length scale. 
The ferromagnetic transition temperature determined within the spin-fluctuation theory is in good 
agreement with experimental values for Gai_ x Mn x As. 



I. INTRODUCTION 

The recent discovery of ferromagnetism in diluted mag- 
netic semiconductors (DMS) such as Gai_ x Mn x Asji*2i&4 
with transition temperature T c ~ 110 K for Mn concen- 
tration x ~ 5%j£ and ~ 150 K in films for x in the range 
6.7 — 8.5%r^ has led to intensive efforts to increase T c 
in view of potential technological applications. Tremen- 
dous interest has also been generated in the novel ferro- 
magnetism exhibited by these systems in which magnetic 
interaction between localized spins is mediated by doped 
carriers. The long-range oscillatory nature of the carrier- 
mediated spin couplings results in a variety of interest- 
ing features of the ferromagnetic state, such as significant 
sensitivity of spin stiffness and transition temperature T c 
on carrier concentration, competing antiferromagnetic in- 
teraction and noncollincar ordering, spin-glass behaviour, 
spin clustering and disorder- induced localization etc., as 
recently reviewed m& 

DMS such as Gai_ x Mn x As are mixed spin-fermion sys- 
tems in which the S = 5/2 Mn ++ impurities replace 
Ga +++ , thereby contributing a hole to the semiconduc- 
tor valence band. However, large compensation due to 
As antisite defects reduces the hole density to nearly 
10% of Mn concentration, which plays a key role in the 
stabilization of long-range ferromagnetic order, and also 
provides a complimentary limit to Kondo systems. The 
interplay between itinerant carriers in a partially filled 
band and the localized moments is conventionally stud- 
ied within a diluted ferromagnetic Kondo lattice model 
(FKLM), wherein — JS/.crj represents the exchange in- 
teraction between the magnetic impurity spin Sj and the 
electron spin crj. 

In order to obtain deeper understanding of the role of 
dilution in terms of the perturbation to carrier spin bands 
due to spatially varying impurity field as well as impurity- 
ficld-induced Zeeman splitting, which control the long- 
range oscillatory nature of the impurity-spin couplings 
through the carrier spin polarization x°(q, u>), we will 
consider a diluted FKLM with ordered impurity arrange- 
ments in this paper. Use of k-space representation allows 
for much larger lattices, thus permitting a more refined 
study of the competing spin interactions with respect 
to carrier concentration, impurity separation, interaction 
strength, and wave vector. More importantly, it provides 



for a quantitative comparison with different theoretical 
pictures, which we briefly review below. 

Long range ferromagnetic interaction between the im- 
purity spins is mediated, in the mean-field (Zcncr model) 
picture ^*Si2iifliiiii^ by a uniform carrier spin polarization, 
which is caused, in turn, by a uniform site-averaged impu- 
rity magnetic field within the virtual crystal approxima- 
tion (VCA). In the weak-field limit (xJS << ep), the car- 
rier spin polarization is proportional to the Pauli suscep- 
tibility xp, and the transition temperature (T c ~ xJ 2 \p) 
is therefore proportional to the Mn concentration x, J 2 , 
the carrier effective mass m* , and N(ep) ~ p 1 ^ 3 , where p 
is the hole concentration. 

An alternative mechanism for the coupling between 
impurity spins involves the carrier-mediated Rudcrman- 
Kittel-Kasuya-Yosida (RKKY) interaction^ The local 
magnetic field Bj = JSj of an impurity spin at site 
J polarizes the electrons locally, and the mobile band 
electrons spread this magnetic polarization in a charac- 
teristic manner: m.j = x?jBj, where x°u represents the 
carrier spin polarization. Another impurity spin Sj cou- 
ples to this local electronic magnetization, resulting in 
an effective RKKY interaction J 2 x5jSj.Sj. Within a 
mean-field (MF) treatment of the resulting Heiscnbcrg 
model^ a similar behaviour for the transition tempera- 
ture (T c ~ xJ 2 xp) is obtained. 

Within the generalized RKKY theory involving the 
non-linear magnetic response^ the carrier-induced spin 
couplings </ 2 X/j(</) go through a maximum with respect 
to both carrier doping concentration p and the fermion- 
spin interaction strength J. This optimization behaviour 
can be qualitatively understood in terms of a competi- 
tion between the increasing magnitude of the carrier spin 
polarization x°ij{J) an d the increasing rapidity of its os- 
cillation, which limits the growth of the spin couplings. 
Similar behaviour was observed in the diluted Hubbard 
modeli^ii^ for the effective spin coupling C/ 2 Xj j{U). 

Long-wavelength magnon excitations provide a com- 
posite measure of the carrier-induced spin couplings in 
the ferromagnetic state, with vansihing spin stiffness sig- 
nalling instability due to competing antiferromagnetic 
(AF) spin interactions. Magnon excitations as func- 
tion of electron density n in the conduction band and 
the spin-fermion coupling J have been studied within 
the concentrated FKLM having a magnetic impurity 
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FIG. 1: An ordered arrangement of impurity spins (•) on a 
square host lattice (o). Also shown are the sublattice labels 
corresponding to the four distinct sites within the unit cell. 



at every lattice site, in the context of heavy fcrmion 
materials^ ferromagnetic metals Gd, Tb, Dy, doped 
EuXil and manganitesiiSiiSiSSiSi Magnon dispersion has 
also been obtained in the context of DMS, where the 
spatially varying impurity field has been treated ex- 
actly numerically^ within the VCA where a uniform 
impurity-induced Zeeman splitting of the carrier spin 
bands is assumed^ and within the coherent potential 
approximation (CPA)M^ Magnon spectrum and tran- 
sition temperature have also been obtained recently for 
Gai_ x Mn x As and Gai_ x Mn x N by means of an effective 
Heisenberg model, whose exchange parameters are ob- 
tained from first-principle calculations^ 

A sublatticc-basis representation for ordered impu- 
rity arrangements is introduced in section II for two- 
dimensional (2d) and three dimensional (3d) cases, and 
the formation of sub-bands due to spatially varying im- 
purity field is discussed. The two energy scales in the 
coupled spin-fermion problem are also introduced, and 
temperature-dependence of impurity magnetization is 
obtained sclf-consistently. Spin-wave excitations are ob- 
tained using the Holstein-Primakoff transformation in 
the large S limit (section III). A detailed study of the 
behaviour of spin stiffness with respect to dimension- 
ality, carrier concentration, interaction strength as well 
as impurity separation (dilution) is presented in sec- 
tion IV. Ferromagnetic transition temperature obtained 
within the renormalized spin-fluctuation theory is com- 
pared with RKKY and VCA results in section V. 



II. DILUTED KONDO LATTICE MODEL 



sity a i = ^l\\cr\^li, where the fermionic field operator 
= ( ) . The host on-site energy is set to zero. 

A. Sublattice-basis representation 

We consider ordered (superlattice) arrangements of 
magnetic impurities on square and cubic host lattices 
with impurity spacings 2a and 3a. Translational symme- 
try in the sublattice basis convenientally allows Fourier 
transformation to momentum space. For concreteness, 
we consider a square host lattice in the following, with 
magnetic impurities placed at every other host site, corre- 
sponding to superlattice spacing 2a and impurity concen- 
tration x = 25%. There are four sublattices, numbered 
a = 1, 2, 3, 4, corresponding to the four sites in the unit 
cell, as shown in Fig. 1. We choose length and energy 
units such that the lattice spacing a = 1 and the hopping 
term t = 1. 

In the MF approximation, the interaction term 
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represents a magnetic coupling of impurity and fermion 
spins with the self-consistcntly determined fermion and 
impurity fields ^(07) and ^-(Sj) = B/, respectively. As- 
suming, without loss of generality, a uniform impurity 
field B/ = Bz 7 Fourier transformation within the sublat- 
tice basis yields the following MF Hamiltonian for spin-cr 
fermion 
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where e^ v) = 2t cos k x r y ) is the hopping energy along the 
x(y) direction and a =+1 and —1 for spin up and down, 
respectively. Here the field operator \&k = (a£ a \ a \ a t) 
defines the sublattice basis, where a£ refer to the fermion 
operator for sublattice index a. The MF Hamiltonian 
is numerically diagonalized to obtain the four eigenval- 
ues EZ w , corresponding to the four sub-band indices 
jU = 1,2,3,4, and the four-component eigenvectors 
yield the fermion amplitudes on the sublattice a. 



We consider a diluted ferromagnetic Kondo lattice 
Hamiltonian 

J , (5.CT I 

on square and cubic lattices with positive hopping t 
between nearest-neighbour (NN) sites i,i + S, which 
yields a host valence band. The second term represents 
the ferromagnetic exchange interaction between local- 
ized impurity spins at site / and the fermion spin den- 



B. Quasiparticle density of states 

The quasiparticle density of states (Fig. 2) shows the 
formation of four and eight sub-bands for the square and 
cubic lattices, respectively, corresponding to the number 
of sublattices. The reflection symmetry about the origin 
of the energy axis reflects the magnetic coupling. The 
sub-band corresponding to impurity sublattice splits off 
and forms an impurity band at energy ~ ±B, which nar- 
rows with increasing B due to decreasing effective hop- 
ping t e ff ~ t 2 1 B between impurity sites. In the undoped 
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FIG. 2: Spin-resolved density of states (DOS), showing the 
formation of sub-bands and impurity bands due to spatially- 
varying impurity field, for (a) d = 2, x = 25% and (b) d = 3, 
x = 12.5%. The impurity field B = JS/2. 



case, both bands are completely filled, and carrier dop- 
ing introduces holes at top of the spin- J, impurity band. 
With further hole doping, Fermi energy moves into the 
spm-| band, particularly for small B. 



C. Impurity magnetization 

With only one energy scale, the NN Heisenberg fer- 
romagnet yields a transition temperature which differs 
from the MF value only by an 0(1) number. However, 
with long-range ferromagnetic order originating from the 
carrier-induced, oscillating impurity-spin couplings, the 
transition temperature for DMS is very different from 
the MF value, which is determined from the bare spin- 
fermion exchange interaction J, and physically repre- 
sents a moment-melting temperature rather than a spin- 
disordering temperature. In order to highlight the order 
of magnitude difference arising from competing interac- 
tions and spin softening, we have also obtained T^ F in 
the MF approximation, which is interesting in its own 
right, and studied recently for the x = 1 case^ 

In the MF approximation, the coupled spin-fermion 
problem is sclf-consistcntly solved from the finite- 
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FIG. 3: Variation of impurity and carrier magnetization with 
temperature, showing a first-order transition for the 3d di- 
luted case (x = 12.5%). Here host bandwidth W = 12t = 10 
eV. 
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of fermion and impurity spins, where /fd is the Fermi- 
Dirac distribution function and (3 = 1/k^T. 

Consider the x = 1 case, for simplicity. There are two 
distinct temperature scales associated with the fermion 
field (erf) ~ \p and the impurity field \ {S\) ~ ^S, re- 
sulting in a characteristic concave (Sf) vs. T behaviour 
in the intermediate temperature regime Jp <C k^T <C 
JS, wherein the nearly free impurity spins yield a typi- 
cal paramagnetic (~ 1/T) response to impurity magne- 
tization. In the diluted case, minority-spin holes prefer- 
entially sit on impurity sites, resulting in an enhanced 
fermion magnetization (erf) ~ p/x, which blurs the dis- 
tinction between the two temperature scales. Conse- 
quently, the concave behaviour in the temperature de- 
pendence of impurity magnetization [Fig. 3] becomes 
prominent only at low doping concentration p and large 
J. With increasing temperature and decreasing impurity 
field, the rapid decrease in (erf) due to redistribution of 
holes from impurity to host sites results in a first-order 
transition, in contrast to the continuous transition ob- 
tained for the x = 1 case. 



III. SPIN- WAVE EXCITATIONS 

Transverse spin fluctuations are gapless, low-energy 
excitations in the broken-symmetry state of magnetic 
systems possessing continuous spin-rotational symmetry. 
Therefore, at low temperature they play an important 
role in diverse macroscopic properties such as existence 
of long-range order, magnitude and temperature depen- 
dence of the order parameter, magnetic transition tern- 
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perature, spin correlations etc. In the following we con- 
sider T = and obtain the spin-wave excitations at the 
non-interacting level. Finite impurity concentration, in- 
teraction strength J, and dynamical effects are treated 
exactly in this approach. 

Applying the approximate Holstcin-Primakoff trans- 
formation in the large S limit 

S+ = biV2S 
Sy = b\V2S 

Sf = S-b\bi (6) 



from the spin- lowering (Sj ) and spin-raising (Sf) oper- 
ators to boson (magnon) creation and annihilation oper- 
ators b\ and bi , which neglects magnon interaction terms 
of order 1/2S, the Kondo-lattice Hamiltonian reduces to 



b\aj 
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where af = <rf ± %a\ . 

In the MF approximation, (&}) = (6/) = (b\bi) = 
and the Hamiltonian decouples into a fermion part 

^fcrmion = a }, a a i+&,<? 7T } ■ a I ( 8 ) 



= 1 E a \,<T ai + s >< 

and a boson part 
nl son = ^(a!)b\bj. 



(9) 



We now obtain the time-ordered, transverse spin- 
fluctuation (magnon) propagator 



xtj-(t-t')=i(* G \T[b I (t)tf J (t')]\* c 



(10) 



within the random phase approximation (RPA), by sum- 
ming over all bubble diagrams 




where the particle-hole bubble in q, uj space 
X°(q,w) 
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FIG. 4: Spin-wave dispersion along main symmetry directions 
of the magnetic Brillouin zone for the (a) 2d (£=25%) and (b) 
3d (a;=12.5%) cases. 



involves integrating out the fermions in the broken- 
symmetry state. It is the particle-hole bubble x°(q, u>) 
which mediates the carrier-induced impurity spin cou- 
plings in the ferromagnetic state, and the oscillatory, 
long-range nature of the spin couplings is effectively con- 
trolled by the Fermi wavevector (hp ~ p 1 ^) and the 
impurity-ficld-induccd Zccman splitting of the carrier 
spin bands. 

The magnon propagator can then be expressed as 



x + (q,w) 
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CJ- (So- £q(w)) 

where the site-diagonal zeroth-order magnon propagator 
1 1 



Ih 



(13) 



cor- 



involves the magnon on-site energy £q 
responding to the energy cost of a spin deviation, and 
£q(w) = =2^x°(q, u>) is the delocalization energy due to 
magnon hopping associated with spin couplings Jjj = 



The energy cost £q of creating a local spin 



deviation is exactly offset by this delocalization-induccd 
energy gain for q, lj — 0, consistent with the Goldstonc 
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FIG. 5: Variation of spin stiffness with (a) hole concentration 
p and (b) exchange energy J for the 2d case (x—25%). 
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FIG. 6: Variation of spin stiffness with hole concentration p 
for the 3d case (x=12.5%). 



mode. Poles in the magnon propagator yield the magnon- 
mode energies by solving 



— Sq — fq(Wq) . 



(14) 



The spin-wave dispersion along the main symmetry di- 
rections is shown in Fig. 4 for the 2d and 3d cases. 



IV. SPIN STIFFNESS 



The spin stiffness D = uj^/q 2 for long-wavelength 
modes provides a composite measure of the impurity spin 



FIG. 7: Comparison of spin stiffness for different impurity 
spin dilutions — (a) x = 1 and 1/4 and (b) x = 1/4 and 
1/9, showing the remarkable enhancement with dilution at 
low doping. 



couplings and stability of the ferromagnetic state. Vari- 
ation of spin stiffness with doping concentration p and 
interaction strength J is shown in Figs. 5 and 6 for 
the 2d and 3d cases. The spin stiffness exhibits opti- 
mization behaviour with respect to both p and J, which 
can be understood within the generalized RKKY theory 
in terms of a competition between increasing magnitude 
of carrier-spin polarization and increasing rapidity of its 
oscillation pi^ An interesting feature of Figs. 5(a) and 6 is 
that due to impurity-host decoupling for higher J values, 
the instability of the ferromagnetic state shifts towards 
" half- filling" (p — ► x) where AF ordering is favoured. For 
a fixed p, a minimum J is required to stabilize the fer- 
romagnetic state, as seen in Fig. 5(b), and the 1/ J- type 
behaviour for large J is associated with narrowing of the 
impurity band due to decreasing effective impurity hop- 
ping t cS - {t 2 /J). 

Fig. 7 highlights the remarkable effect of dilution 
on spin stiffness — dilution actually enhances the spin 
stiffness at low hole doping, demonstrating an optimiza- 
tion with respect to impurity concentration as well. The 
above result highlights the key requirements of impurity 
spin dilution as well as low hole doping concentration 
for a robust ferromagnetic state, both of which are met 
in the heavily compensated material Gai_ x Mn x As, with 
x ~ 5%. It is the reduced particle-hole energy gap in the 
fermion polarization bubble x°(q, w) for the diluted case 
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FIG. 8: Comparison of spin stiffness for the 3d case with 
(x — 1/8) and without (x = 1) dilution for same J, and with 
the VCA result (x = 1, J = 10/8). 
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FIG. 9: Comparison of the transition temperature T c ob- 
tained from the renormalized spin-fluctuation theory (SF) 
with the RKKY and MF results for the 3d case (x=l2.5%). 



which relatively enhances the oscillating carrier-spin po- 
larization, resulting in greater spin stiffness for low dop- 
ing and higher sensitivity to competing interactions with 
increasing doping. Indeed, the limit of vanishing energy 
gap (RKKY theory) yields an upper bound to the spin 
stiffness and transition temperature (see Fig. 9). 

Fig. 8 shows the enhancement of spin stiffness with 
dilution at low doping for the 3d case. Also shown is a 
comparison of the spin stiffness with the VCA result for a 
uniform impurity field B = xJS/2 with identical lattice- 
average. The VCA result provides a fair approximation 
in the low doping regime, which can be understood in 
terms of an averaged impurity field seen by doped carriers 
with wavelength longer than the impurity spacing. 



V. TRANSITION TEMPERATURE 

Determination of the magnon spectrum in the carrier- 
induced ferromagnetic state allows for an estimation of 
the Curie temperature T c in three dimensions. As the 
ferromagnetic state is characterized by small spin stiff- 
ness due to competing interactions, the dominant contri- 
bution to reduction in magnetization is from the thermal 
excitation of long-wavelength magnon modes. Therefore, 
T c can be estimated in terms of an equivalent Heisenberg 
model with matching spin stiffness. For a spin-iS Heisen- 
berg ferromagnet with nearest-neighbour interaction J 
on a simple cubic lattice (coordination number z~6), the 
magnon energy 

w q = zJS{\ - 7 q ), (15) 

where 7 q =(cos2q , ;c -|-cos2(jfj,-|-cos2g z )/3, corresponding to 
the magnetic lattice spacing 2. Indeed, the magnon en- 
ergy is maximum for q x = q y = q, = it/ 2, as also seen in 
Fig. 4 (at R). Considering the small q limit, we obtain 
the spin stiffness D = uj^/q 2 = 4JS = 10J for S = 5/2. 
Within the renormalized spin-fluctuation theory, the 



transition temperature 

T? F - ^^ /sp 1 = \d{S + l)/ s -; (16) 

where the spin-fluctuation factor /sf = (1 — 7 q ) _1 ~ 
1.5 for the cubic lattice. As T C SF < T C MF , the impu- 
rity mean field is essentially unchanged at low tempera- 
ture; therefore we have assumed the zero-temperature 
value for the spin stiffness. In the RKKY- VCA ap- 
proach, the lattice-averaged impurity field is assumed 
to decrease with impurity magnetization and vanishes 
as T — > T c iSiS4 Considering a bandwidth of the or- 
der of that of GaAs (W = 12t - 10 cV), we obtain 
T C SF = (35/36)L> in eV. The calculated T C SF (Fig. 9) 
is in qualitative agreement with experimental results for 
Gai^MnzAs. 

Figure 9 also provides a comparison with the RKKY 
result^ 

T RKK Y = ;rj2 5(5+1)^ 

12fc B 

obtained within the mean-field approximation from the 
carrier-induced RKKY spin couplings Ju = J 2 x°ij- I n 
the limit of low doping concentration (p <^ 1), the host 
density of states 

N(E F ) = ^(3^ 2 P) 1/3 , (18) 

yielding the characteristic p 1 / 3 behaviour of T^ KKY . It 
is seen in Fig. 9 that the RKKY result drastically over- 
estimates T c , mainly due to neglect of the impurity- 
ficld-induccd energy gap in the fermion polarization bub- 
ble, and also fails to capture the optimization behaviour 
with carrier doping, impurity spin dilution, or interaction 
strength. The MF transition temperature T^ F , which 
physically represents the moment-melting temperature, 
is an order-of-magnitudc still higher. 
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VI. CONCLUSIONS 

Magnon excitations in the carrier-induced ferromag- 
netic state of the diluted FKLM were studied for sev- 
eral ordered impurity arrangements. Use of momentum- 
space representation allowed for a refined study with re- 
spect to variation of doping concentration, interaction 
strength, impurity spin dilution, and wave vector. The 
spin stiffness was found to exhibit a characteristic op- 
timization behaviour with respect to hole doping con- 
centration, interaction strength, as well as impurity spin 
dilution, which can be qualitatively understood in terms 
of a competition between increasing magnitude of the 
carrier-spin polarization and increasing rapidity of its os- 
cillation. 

The exact treatment of impurity spin dilution and 
interaction strength allowed for a detailed comparison 
with results of existing theories such as MF, RKKY, and 
VCA. We found that the RKKY result drastically over- 



estimates T c due to neglect of the impurity-ficld-induced 
Zeeman splitting of the carrier spin bands, and also does 
not capture the optimization behaviour with carrier dop- 
ing, impurity spin dilution, or interaction strength. The 
VCA approach involving a uniform lattice-averaged im- 
purity field was found to provide a fair approximation for 
spin stiffness in the low doping regime, but fails to de- 
scribe the competing interactions at higher doping when 
the hole wavelength becomes comparable to impurity 
spacing. 

The enhancement of spin stiffness with dilution at low 
doping concentration due to the reduced impurity-field- 
induced energy gap in the fcrmion polarization bubble 
highlights the key requirements of impurity spin dilu- 
tion as well as low hole doping concentration for a robust 
ferromagnetic state of diluted magnetic semiconductors, 
both of which are realized in the heavily compensated 
material Gai_ x Mn x As, with x ~ 5%. 
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